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Abstract 

Charged scalar particles introduced in some extensions of the standard model can 
induce r leptonic decay at tree level. We find that with some charged SU(2)-singlet 
scalar particles, like ones introduced in Zee-type models, r leptonic decay width is 
always smaller than what is predicted by the standard model, therefore they may offer 
a natural solution to r decay puzzle. To be more specific, we examine some Zee-type 
models in detail to see if at the same time they are acceptable in particle physics, 
cosmology and astrophysics. It is shown that r decay data do put some constrains on 
these models. 



PACS numbers: 13.35.+S, 13.10. +q, 12.15.Ff 



r lepton is an interesting system to test the standard model (SM) and search for new 
physics, since r is the heaviest lepton yet known. There is a long-standing puzzle in r 
lepton decays, which has received attention for some years. The puzzle is that the measured 
r lifetime may be longer than one expected in the SM with three families |2[. From 
Particle Data Group (PDG) [|, the measured r lifetime is r exp = (3.05 ± 0.06) x 10~ 13 s, 

+2.7 

while the SM's expectation is r th = (2.87 ± 0.07) x 10~ 13 s, where m T = 1784.1 MeV 

-3.6 

is used. So one sees that the measured r lifetime is about 2.3<r higher than SM expectation 
value ||. The latest measurement of r mass at BES Q m T = 1776.9 ±0.5 MeV somehow 
relaxes the r lifetime problem. But this downward shift of m T is not enough. r exp is still 
about 1.9a higher than r th (= 2.92 ± 0.04 x 10~ 13 ) . Of course it is very possible that this 
r decay puzzle will disappear when new measurements of leptonic decays become available, 
as the expected value does not deviate too much from the measured one (In fact, it is 
noticed that there are some new measurements after PDG || on r lifetime and leptonic 
decay branching ratios. We will comment on that at the end of the paper). However, we 
feel that there are some theoretical motivations for taking this puzzle seriously, such as the 
existence of a fourth generation or charged scalar particles (the latter case will be discussed 
in detail later). 

One simple solution to r decay puzzle is to introduce a fourth generation |2|, |5|]. If the 
mixing between r neutrino v T and the fourth heavy neutrino (it must be heavier than 45.3 
GeV from LEP Z-width constraints ||) is around sin 2 9 mix ~ 0.05, the central value of 
T exp j g j n consistency with the corresponding theoretical expectation value in this model. 
We denote g e , and g T as weak couplings of e, fi and r leptons respectively, universality 
of weak interaction means g e = g M = g T . Any deviation from this relation (for example 
9e — 9n 7^ 9t) implies a violation of the universality (of r). In the four generation model, 
r universality in the charged and neutral current sector is obviously violated by a small 
amount 0. Nevertheless, this is not favored by the neutral-current data from Z decay 
which agrees with the universality of the weak interaction of e, [i and r leptons at the 
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level of precision better than 0.5% ||. Another simple solution assuming a mixing of v T 
with a singlet neutrino v || is also in conflict with this Z decay data. In addition, following 
Marciano's analysis || but using m T = 1776.9±0.4±0.3 MeV, we see that for r semileptonic 
decays r — > v T ir/K and r — ► z/ r 7r _ 7r experimental values agree with the SM's prediction 
very well (within la). So this is another evidence supporting r universality in semileptonic 
decay. 

In this letter we discuss the effects induced by scalar particles in r leptonic decay. With- 
out losing generality we will consider SU(2)-singlet, doublet and triplet scalars. Generally, 
for a SU(2)-singlet or triplet scalar h a b, the couplings between scalars and leptons may be 
introduced as AL y = f a blJbh a b, for SU(2)-doublet scalars AL y = f a J-a^b^ab, where l a is lep- 
ton doublet, e£ is lepton singlet and a, b denote family indices. Due to the fermi statistics, 
fab is antisymmetric and symmetric for singlet and triplet scalars respectively. While for 
doublet scalar there is no any constraint on the structure of f a b- With these interactions, it 
is easy to see that for r leptonic decay, for example r — ► iw T v^ two processes contribute at 
tree level. One is SM's W-boson exchanging process, the other is through exchange of the 
h a b particle (Fig. 1). The general features we find are that for singlet h a b, the interference 
term between these two processes are always negative; for triplet h a b, if a ^ b, the interfer- 
ence term is always positive, if a = b and f aa fbb > the interference term is negative. These 
are because of the antisymmetric and symmetric properties of f a b- As for the doublet, if we 
take all of the fermion masses in final states as zero, then the interference effect vanishes. 
This is a result of the fact that the Yukawa coupling changes the chirality of the leptons, 
and therefore the interference term must be proportional to the charged lepton mass in the 
final states. Consequently, we see that with the singlet or triplet scalar couplings (with 
a = b and faafbb > 0) r leptonic widths are smaller than that predicted in the SM, with 
triplet scalar couplings (in the case of a ^ b) or doublet scalar couplings r leptonic decay 
widths are larger than that in the SM. Hence r decay data suggests that one should consider 
models with singlet or triplet (with a = b and f aa fbb > 0) scalars. We show that the scalar 
particle introduced in Zee-model has exactly the required property . So as a good example, 
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we will discuss r decay puzzle concentrating only on Zee-type models. 



Zee model was proposed to generate Majorana neutrino masses [ID|. Recently, it was 
found that it can generate large neutrino transitional magnetic moment, while keeping 
neutrino masses small naturally pl| . Zee- type models are also proposed to incorporate 



neutrinos with mass at the order of KeV |T2| and at the same time to give a solution to 



solar neutrino problem (SNP) ||14|| . In this sort of models, some charged scalar particles h are 
introduced, which carries two units of lepton charge. The point is that the r leptonic decay 
widths in these models are always smaller than that predicted in SM. Also r universality 
in the neutral current sector is not violated. On the other hand, since h does not couple to 
quarks, we expect the r universality is well respected in semileptonic decays. These are in 
perfect agreement with the Z decay and r leptonic decay data as well as the more precisely 
measured value Ig^/gel = 1.0031 ± 0.0023 obtained from 7r decay 

h a b(a 7^ b) couple to leptons as 

AL y = l/2f ab l T a Cir 2 l b h ab + h.c. (1) 

where C is the Dirac charge conjugation matrix and f ab is antisymmetric due to fermi 
statistics. h ab carries two units of lepton numbers (L a ,L b ). One sees that this interaction 
has a global U{l) e x t/(l) M X U (1) T symmetry. So lepton numbers L a are not violated through 
this interaction. There are three independent couplings f efl , f er and / Mr . We assume that f efJi 
is considerably smaller than f eT and / MT , so that we don't need to readjust fermi constant 
( in fact one or two order of magnitude smaller is enough, since we don't want to fine-tune 
f efJi either). This is also consistent with the constraint set by universality between beta and 
//-decay |TD|, [I5|]. The leptonic decay width (including electroweak radiative corrections of 



the SM) in Zee-type models reads 



G l m r ,u m2 iM^ | 3m rwT | «(m T ) ^25 2 

f2 /i nnn„y/ \2 



C h = 3.0ff T (100GeV/m h y (2) 
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here f(x) = 1 — 8x + 8x 3 — x 4 — 12x 2 \nx, is fermi coupling constant and a is fine 
structure constant. If we consider the weak scale as the only physical scale in this work, 
i.e. rrih — 100 GeV, and take f) T = 0.12, Ti will be about 4% smaller than what predicted 
in SM. Consequently r lifetime is about 4% longer than SM's prediction as the central 
value of r exp implies. It is interesting that with a reasonable choice of coupling constant 
fi T ~ 0.1, t lifetime is about a few percent longer than the SM's prediction. This is a 
general prediction of Zee-type models. Since x is proportional to ff Tl a reduction of fi T 
by one order of magnitude will decrease x by two order of magnitude. In other words, fir 
can be determined very well from r lifetime measurement, if r lifetime is indeed different 
from the SM's prediction. A few percent deviation between experimental measurement and 
the SM's expectation implies fi T ~ 0.1 and an inconsistency at 10 -3 level corresponds to 
fir ~ 0.02. In addition, f^ is directly related to neutrino masses, mixing and magnetic 
moments. Precision measurement of r lifetime will constrain these parameters and then 
affect the descriptions of some phenomena of cosmology and astrophysics. In order to see 
how well the models works when the parameters are fixed from the r decay, we examine 
some of the aspects which are sensitive to the values of f a b in some concrete examples. 



BFZ model |TT| In this model both neutrino masses and magnetic moments are 
generated at two loop level, but the neutrino masses are suppressed by a factor proportional 
to the mass square of the charged leptons due to the spin suppression mechanism. Individual 
lepton number is certainly violated because of the non-zero neutrino transitional magnetic 



moments and masses. This violation is only due to the scalar interaction [17| 



Kft&fafi - r P €) (3) 

here <p~, (j)° a belong to Higgs doublets. The neutrino mass matrix for three lepton flavors 
reads 

(n _ \ 

(4) 
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and m a b oc fabi^l — m 1)-> so one nas m ^ = m << M = m eT ~ m^ T . Therefore this mass 
matrix has an approximate symmetry L e + — L T . The eigenvalues of this mass matrix 
are m\ ~ m, 7712,3 ~ V2M ± m/2 indicating the mixing angle between v e {y^) and v T 
to be order of m/M, and u e , are mixed with large mixing angle close to 45°. Taking 
fer ~ ftir ~ 0.1 and / eA1 < 10~ 2 to satisfy f efM « f er , / MT , one has neutrino transitional 
magnetic moments (d v ) eT ~ {d v )^ T ~ 10 _11 /ib and (d u ) e/1 < 10~ 12 fiB- With m smaller than 
10~ 4 eV, and M as large as 0.1 eV, it is easy to see that the relevant squared mass difference 
(Am 2 ) eM ~ M 2 ~ 10 -2 eV 2 in v e and oscillation is too large and (d u ) e ^ < 10~ 12 fiB is too 
small so that it can not provide a solution to SNP either through neutrino oscillation or 
spin-flavor precession between v e and v^. However, we notice that the oscillation between v e 
and Up can be responsible for the recently reported deficiency of atmospheric (ANP) [ |18|| . 
Because the squared mass difference and mixing angle perfectly fit the required parameter 



range ||18|| . As for the SNP we have a solution resorting to non- resonant spin-flavor precession 
between v e and v T when v e goes through magnetic field inside the Sun. The reason is that 
the magnetic moment (d v ) eT can be sufficiently large ~ 10^ 11 ub and (Am 2 ) eT ~ Mm < 10 -5 
eV 2 is as small as required |19[]. Of course we should also check on all of the experimental 
constraints from particle physics, cosmology and astrophysics. As there are only three light 
neutrinos in our case, it seems that there is no problem with the limits from cosmology 
and astrophysics. Neutrino oscillation experiment gives some restrictions on mixing angle 
and mass difference between two neutrino species. Because of the approximate symmetry 
L e + Lp — L T the only possible disagreement with experimental constraints could happen in 
v e and z/ M oscillation. Given a large mixing angle, neutrino oscillation experiment requires 
(Am 2 ) eM < 0.09 eV 2 0, so our results agree with this restriction. However, if we take a 
much more stringent limit (Am 2 ) e(i < 1.5 x 10~ 3 eV 2 [2D|, then our prediction is in conflict 
with this limit. Surely we have some freedom to tune the parameters to make M a few 
times smaller in order to satisfy this limit. But at the same time we also reduce magnetic 
moment (d u ) etl by a same factor. This is not what we would like to do. Another stringent 
constraint comes from the measurement of rare \i decay /i — > e^y 0. In present case this 
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decay can happen but is dominated by two loop diagrams plR , as lepton number is violated 
only through scalar interaction (3). A rather conservative estimate gives / 2 < 10~ 2 , this is 
consistent with our requirement f efl << 10 _1 . 



BH model Jl2j This model is based on a global lepton flavor symmetry G = U(l) e x 
U{\)n x U(1) T and G is spontaneously broken down to £/(l) e _ M+r at weak scale. Because of 
the U(l) e -n +T symmetry, neutrino mass matrix has the form 

fn __ „ \ 

(5) 
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m efl and m TfJ , arise from one loop diagrams [14], giving M = m e[ii ~ l/167r 2 / MT A T m T , m = 
m eT ~ l/167r 2 / eAt A e m M , where A's are scalar coupling constants and we expect these to be 
of the same order of magnitude. Solving the eigenvalues of this matrix, one gets a massless 
neutrino which is mostly v e mixed with v T by a mixing angle Q$ — m/M, and a massive 
ZKM neutrino with mass ~ M. The latest measurement on searching for 17 KeV neutrino 
sets a restriction 9s < 1CT 3 with 95% CL [p3|| . This limit requires f ef i/f^T < 1CT 2 , which 
is again consistent with our assumption on f e/I . The mass of the heavy ZKM neutrino is 
naturally about 10 ~ 100 KeV. To avoid any trouble in cosmology and astrophysics , this 
heavy neutrino must decay sufficiently fast and the dominant decay is through v T — > u e F, 
where F denotes flavons associated with the spontaneous breaking of G. The interesting 
prediction of this model is also on r decay, i.e. Bir — > eF) ~ 10 -4 . With our choice of 
parameters, we get same value for this branching ratio. If we demand the lifetime of the 
heavy neutrino is less than ~ 10 5 sec to be consistent with the conventional mechanism 
for large scale structure formation [12], then we may put a lower limit 10~ 4 ~ 10~ 5 on the 
value of f efl . So in this model we can more or less fix the parameters f a b- As the global 
symmetry G is spontaneously broken, there may be some additional tree level r leptonic 
decay processes, like r — > /iz/^z/g, but this kind of processes are always proportional to the 
f eil . Therefore the corresponding branching ratios are negligibly small. 
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Extension of BH model In order to incorporate SNP and ANP, BH model was 
extended to four neutrino species including one sterile neutrino n [O, [IJJ. The symmetry 



group is extended to G — U{l) e x £/(l) M x U(1) T x U(l) n . The four neutrino species may 
make up one Dirac neutrino and one ZKM neutrino. Generally there are two possibilities to 
incorporate SNP. One is that v e and are combined to form a light ZKM neutrino, v T and 
n (or n c ) make up a heavy Dirac neutrino with mass around 10 KeV. Applying the BFZ 
mechanism on the light ZKM neutrino part, one can generate a large transitional magnetic 



moment between v e and u^, hence providing a spin-flavor precession solution to SNP |{U 
However, in present case BFZ mechanism does not work since f efl is very small, and therefore 
the expected neutrino transitional magnetic moment is not large enough. Another possibility 



is that and v T are combined to form a heavy ZKM neutrino, and v e and n (or n c ) make up 
a light Dirac neutrino. Both SNP and ANP can be solved naturally in this scheme resorting 
to the Planck scale effects 0, |16| . As in BH model the mixing angle 9 s between v e and the 
heavy neutrino could be naturally smaller than 10 -3 due to the smallness of f efl . In addition, 
with 9s < 10~ 3 the decoupling temperature of the sterile neutrino n in the early universe 
is higher than the QCD phase transition temperature, consequently the effective number 
of neutrino species at the time of nucleosynthesis is smaller than 3.3 |TJ|]. Also, there are 



more r leptonic decay channels in this case than in BH model, for example r — > ev^n or 
t — > ePfj,(n ) could happen at tree level. But the amplitudes of the tree diagrams for these 
channels are proportional to light neutrino mass, so they are too small to give rise to any 
observable effect. 

In summary, it is shown from above discussion that in Zee-type models the branching 
ratios of the r leptonic decays are naturally smaller than the predictions of the SM. Therefore 
t lifetime in these models is longer than that in the SM. This suggests a solution to r decay 
puzzle. Moreover, the universality between e, /x and r is not violated in the neutral current 
sector and is also respected in semi-leptonic decays of r. These are favored by current 
experimental data. 
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We discussed three Zee-type models in detail, but only concentrated on the issues which 
could be sensitive to the values of f a b fixed by solving r decay puzzle. As to other issues, 
there are lengthy discussions in original papers on these models. We have analyzed three 
specific examples to illustrate the idea of how to incorporate r decay puzzle and take into 
account other related issues of particle physics, cosmology and astrophysics. It is very 
possible that there are some other kind of models which can do the same job |[22|| . We used 
Zee-type models as an example, since we think that the models discussed above are among 
the most popular and simplest extensions of the SM . 

All of the previous analysis are based on the r decay data which implies existence of 
the r decay puzzle. However, whether there is really a disagreement between experiment 
and the SM on r lifetime is not very clear yet. In this work we used the data from PDG 
0, however, meanwhile some new data have become available. In ref. ||, it is reported 
that the new world average data of r decay agree with the SM within la. Nevertheless the 
most precise measurement among these new experiments on r leptonic decay from CLEO 
indicates the discrepancy is still around 2a [24, 25]. 

Anyway, further efforts on precision measurement of the r decay is certainly very much 
desirable. A confirmation of the r decay puzzle will imply some new physics beyond the SM, 
probably as suggested by the models with some singlet charged scalar particles, like Zee-type 
models, or models with triplet scalars. A negative result is of course another evidence in 
supporting the SM and will furthermore constrain the parameter space of these models. 
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Figure Caption 



Fig.l Tree level diagrams which contribute to r leptonic decay with additional 

scalar particles. 
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